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Abstract 



We reconsider the minimal supersymmetric SO(10) model, where only one 10 and one 126 
Higgs multiplets have Yukawa couplings with matter multiplets. The model is generalized 
to include CP-violating phases, and examined how well its predictions can meet the 
current neutrino oscillation data. Using the electroweak scale data about six quark masses, 
three angles and one CP-phase in the Cabibbo-Kobayashi-Maskawa matrix and three 
charged-lepton masses and given tan (3 (the ratio of vacuum expectation values of a pair of 
Higgs doublets), we obtain the Pontecorvo-Maki-Nakagawa-Sakata matrix and the ratio, 
Amg/Am^, as functions of only one free parameter in the model. In our analysis, one- 
loop renormalization group equations for the gauge couplings, the Yukawa couplings and 
the effective dimension-five operator are used to connect the data between the electroweak 
scale and the grand unification scale. Fixing the free parameter appropriately, we find, 
for example, sin 2 26 12 ~ 0.72, sin 2 2 6 23 ~ 0.90, sin 2 26 13 ~ 0.16 and Ami/Am^ ~ 0.19 
with tan/3 = 45, which are in agreement with the current neutrino oscillation data. 
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1 Introduction 



Supersymmetry (SUSY) extension is one of the most promising ways to provide a solution 
to the gauge hierarchy problem in the standard model 0. The minimal version of this 
extension is called the minimal supersymmetric standard model (MSSM). Interestingly, 
the experimental data support the unification of the three gauge couplings at the scale 
M G ut ~2x 10 16 GeV with the MSSM particle contents j2J. At high energies, our world 
may be described by a SUSY grand unified theory (GUT) with a simple gauge group such 
as SU(5) or SO (10) into which all the gauge groups in the standard model are embedded 
and unified. 

On the other hand, much information about quark and lepton mass matrices has been 
accumulated in these decades. Especially, there exist, at present, strong evidences of neu- 
trino masses and mixings through the interpretation of the (active) neutrino oscillations 
as the solutions to the solar neutrino deficit [HI IU El El III E] and the atmospheric neutrino 
anomaly 9 r According to the results, the MSSM has to be extended so as to incorporate 
neutrino masses and mixings. 

In these points of view, the supersymmetric SO(10) GUT model is one of the most 
well-motivated theories, since the right-handed neutrinos are naturally incorporated and 
unified into the same multiplet, 16, together with the ordinary matters in the MSSM. 
Furthermore, this model has an additional advantage that the smallness of the neutrino 
masses can be naturally explained through the see-saw mechanism ^U] with large right- 
handed Majorana neutrino masses. 

There are many possibilities for the Higgs sector that will make the model realistic: 
namely, it correctly breaks the SO(10) gauge symmetry into the standard model one, re- 
produces observed fermion masses, mixings and CP-phase, and realizes the doublet-triplet 
Higgs mass splitting (in terminology of SU(5) GUT) etc. In this paper, we concentrate 
our discussion on the fermion masses and mixings, and consider the so-called "minimal" 
SO(10) model where only one 10 and one 126 Higgs multiplets have Yukawa couplings 
with quarks and leptons. It is very interesting to investigate this model, since it can fix 
all the fermion mass matrices including, at present, ambiguous and undetermined ones. 

The minimal model was first seriously considered by Babu and Mohapatra JT] (non- 
SUSY case) and Lee and Mohapatra J2| (SUSY case), where the neutrino oscillation 
parameters were predicted with inputs of six quark masses, three mixing angles in the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix and three charged-lepton masses. Unfortu- 
nately, the predictions were revealed to contradict the current experimental data. For 
example, the predicted mixing angle, sin 2 2023, in the Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) matrix was found to be too small to be consistent with the the atmospheric 
neutrino oscillation data 9j. It seems to be inevitable that the minimal model should be 
extended by introducing new Higgs multiplets so as to incorporate the realistic neutrino 
oscillation parameters. 

However, note that, in Refs. [TT] and [T2], only the CP invariant case was analyzed. 
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It was recently pointed out ^3] that, once CP-violating phases are taken into account 
in the minimal model, there exits the parameter region which cannot be excluded by 
the current neutrino oscillation data. This observation brings a new hope to obtain the 
minimal SO (10) model compatible with all the observed data of fermion mass matrices. 
Nevertheless, more elaborate studies are needed, since, in jTH], the GUT relation among 
the charged fermion mass matrices was applied to the mass matrices at the electroweak 
scale. This treatment is incomplete because the GUT relation is valid only at the GUT 
scale. In fact, we can see that Yukawa couplings evolve according to the renormalization 
group equations (RGEs), and the GUT relation is lost after the running in the "long 
desert" between the GUT scale and the electroweak scale. 

In this paper, we pursue this program in more correct and complete way. Our strategy 
is the following. First, we evaluate the data of charged fermion mass matrices at the GUT 
scale by extrapolating the data at the electroweak scale according to their Yukawa coupling 
REGs. Next, these GUT scale data are applied to the GUT relation that is generalized to 
include CP-violating phases. This leads us to the explicit form of neutrino mass matrix 
at the GUT scale via the GUT relation. Lastly, running it back to the electroweak scale 
according to the RGE for the effective dimension-five operator , we compare our result 
with the current neutrino oscillation data. 

This paper is organized as follows: in the next section, we give our basic framework of 
the minimal SUSY SO(10) model with one 10 and one 126 Higgs multiplets, and discuss 
the GUT relation among the fermion mass matrices. In Sec. 3, RGEs we use in our 
analysis are introduced. In Sec. 4, numerical analysis is performed, and the results are 
presented. Other predictions from our resultant mass matrices are discussed in Sec. 5. 
The last section is devoted to summary and comments. 

2 Minimal SO (10) model and fermion mass matrices 

We consider the minimal SUSY SO(10) model with a pair of 10+126 Higgs multiplets, 
only which have Yukawa couplings (superpotential) with the matter multiplets such as 

w Y = YSie^oie,- + r ] l 6 i6^ 126 i6 j , (i) 

where 16j is the matter multiplet of the i-th generation, H w and H^q are the Higgs 
multiplet of 10 and 126 representations under SO(10), respectively. Note that, by virtue 
of the gauge symmetry, the Yukawa couplings, Y 1Q and Yi 26 , are complex symmetric 3x3 
matrices. In our model, these Yukawa couplings are assumed to be real as in the original 
model ^1] to keep the number of free parameters minimum. 

In the SUSY GUT scenario, the GUT gauge symmetry is broken at the GUT scale, 
Mqut ~ 10 16 GeV, into the standard model one. There are some different breaking 
chains suitable for our aim. Here, suppose the intermediate pass through the Pati-Salam 
subgroup [T3|, G 42 2 = SU(A) C x SU(2) L x SU(2) R , for simplicity. Under this gauge 
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symmetry, the above Higgs multiplets are decomposed as 10 — » (6, 1, 1) + (1, 2, 2) and 
126 -> (6, 1, 1) + (10, 3, 1) + (10, 1, 3) + (15, 2, 2), while 16 -> (4, 2, 1) + (4, 1, 2) for the 
matter multiplet. Breaking down to the standard model gauge group, SU (4) c x SU (2)r — > 
SU(3) C x U(1)y, is accomplished by vacuum expectation value (VEV) of the (10,1,3) 
Higgs multiplet. Note that the Majorana masses for the right-handed neutrinos are also 
generated by this VEV through the Yukawa coupling Y U q in Eq. (JTJ). 

After this symmetry breaking, we find two pair of Higgs doublets which are in the 
same representation as the pair in the MSSM. One pair comes from (1,2,2) C 10 and 
the other comes from (15, 2, 2) C 126. Using these two pairs of the Higgs doublets, the 
Yukawa couplings of Eq. ((T| are rewritten as 

+ v i (^10-^10 — 31^26-^126) $-3 + e * (^10-^10 — ^li^X2&) 

+ ViiY&v^Vj, (2) 

where u, d, V and e are the right-handed SU{2) l singlet quark and lepton superfields, q and 
£ are the left-handed SU(2)l doublet quark and lepton superfields, H™q and H^q are up- 
type and down- type Higgs doublet superfields originated from H w and Hue, respectively, 
and the last term is the Majorana mass term of the right-handed neutrinos with VEV of 
the (10, 1, 3) Higgs, Vr. Quarks and leptons acquire Dirac masses through VEVs of these 
Higgs doublets. Note that, in that case, the entry of the Clebsch-Gordan coefficient, —3, 
in the lepton sector plays the crucial role so that the unwanted GUT relations, ma = m e 
and m s = m^, are corrected ^T] in the same manner discussed by Georgi and Jarlskog 

Here, remember that the gauge coupling unification succeeds with only the MSSM 
particle contents. This means that only one pair of Higgs doublets remains light and the 
others should be heavy (> Mq). It is highly non-trivial to accomplish this task naturally, 
while keeping the crucial role of H^ to the fermion masses. This problem is essentially 
the same problem as the doublet-triplet Higgs mass splitting problem, and is referred as 
the "doublet-doublet splitting problem". In this paper, we assume that this splitting is 
realized. 

Now let us see the low energy superpotential of our model. After realizing the doublet- 
doublet Higgs mass splitting, we obtain only one pair of light Higgs doublets (H u and Hd) 
in the MSSM, which are, in general, admixture of all Higgs doublets having the same 
quantum numbers in the original model such as 

H u = a u H^ + P U H U Q + • • • 

Hd = OLdH^ Q + f3dHf 2& + ■ ■ • , (3) 

where "+ • ■ •" denotes the mixing with other Higgs doublets in the model (about which 
we do not discuss), a Uj d and (3 u ,d denote elements of the unitary matrix which rotate the 
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flavor basis in the original model into the (SUSY) mass eigenstates. Omitting the heavy 
Higgs mass eigenstates, the low energy superpotential is described by only the light Higgs 
doublets H u and H d such that 

W Y = u t (a u Y^ + p u Y!i 6 )H u q,+d t (a d Y^ + p d Y;i 6 )H dqj 
+ 77, (a u YH - W u YlL) HJi + ei (a d YH - 3(3 d Y^) H d £j 
+ ^(yjLv^Vj, (4) 

where the formulas of the inverse unitary transformation of Eq. (JHJ), #10 = & u,d H Utd + • • • 
and if "26 = P u,d H u ,d + •••■, have been used. Note that the elements of the unitary matrix, 
a u,d anc j pu,d^ are - n g enera } complex parameters, through which CP-violating phases are 
introduced into the fermion mass matrices. 

Providing the Higgs VEVs, H u = v sin/3 and H d = v cos ft with v = 174GeV, the 
quark and lepton mass matrices can be read off as 



M u 


= C10M10 + Ci 26 Mi26 


M d 


= M 10 + M 126 


M D 


= C10M10 - 3ci2 6 Mi2 6 


M e 


= M 10 -3M 126 


M R 


= c R M 126 , 



where the mass matrices are defined as M 10 = Y w a v cos (3 and M 12 6 = Y^ft v cos /3, 
respectively, c w = (a u /a d ) tan/3, ci 2 6 = (P u /P d ) tan/3, and cr = Vn/(f3 d v cos/?). All the 
quark and lepton mass matrices are characterized by only two basic mass matrices, M10 
and M126, and three coefficients. This means the strong predictability of the model to the 
fermion mass matrices in the MSSM. 

Now let us count how many free parameters are included in the above mass matrices, 
not counting cr. Remember that we began with real symmetric matrices Y w and Yi26- By 
flavor rotation of the GUT matter multiplets, we can begin with the basis where M10 is real 
and diagonal, which includes three free parameters. There are seven parameters in M\2§: 
six parameters originate from the symmetric real matrix Y126, arid one from the phase of 
f3 d . Adding the four parameters from complex cio and c\2q, we find 14 free parameters 
in total, which are evaluated as six quark masses, three angles and one CP-phase in the 
CKM matrix and three charged-lepton masses, but one parameter is left free ^Zj. Now 
two free parameters are newly introduced in addition to 12 free parameters analyzed in 
[TT] and ^2], by which our model is generalized to include CP- violating phases. 

The parameter, cr, appears in the right-handed Majorana neutrino mass matrix, by 
which masses of the right-handed Majorana neutrinos are determined after M 12 6 are fixed. 
Since Vr ~ M G , cr ~ vr/v ~ 10 14 GeV is naturally expected. However, we treat cr as 
a free parameter in the following analysis, since the Higgs multiplet breaking G422, in 
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general, mixes with other Higgs multiplets being in the same representation under G^2- 
Also there should exist Higgs multiplets which also have G422 breaking VEV to cancel 
the D-term VEV and ensure SUSY. In addition, threshold corrections would introduce 
an uncertainty of vr. 

The mass matrix formulas in Eq.© leads to the GUT relation among the quark and 
lepton mass matrices, 

M e = c d (M d + kM u ) , (6) 

where 

3Cio + C126 

Cd = 

ClO — c 126 

« = -0 " • ( 7 ) 

3Cio + C126 

Without loss of generality, we can begin with the basis where M u is real and diagonal, 
M u = D u . Since M d is the symmetric matrix, it can be described as M d = V^ M D d V^ M by 
using the CKM matrix Vkm and the real diagonal mass matrix D d . Considering the basis- 
independent quantities, tr[M|M e ], tr[(M ( jM e ) 2 ] and det[M^M e ], and eliminating \c d \, we 
obtain two independent equations such that 

tr[(M e f M e ) 2 ] 

~~ — — 2 ' 

det [MjM e ] 

99 95 ) 
m£ m~ 

where M e = V^ M D d Vj <M + kD u . With the input data of six quark masses, three angles 
and one CP-phase in the CKM matrix and three charged-lepton masses, we can determine 
k and \c d \, but one parameter, the phase of c d , is left undetermined [T7j . Now the mass 
matrices, M w and M 12 6, are described by 

M10 = + (10) 

M 126 = ^wC-^, (11) 

as the functions of a, the phase of c d , with the solutions, |c<f| and k, determined by the 
GUT relation. 

Note that the GUT relation of Eq. © is valid only at the GUT scale. Therefore, 
the input data needed to determine k and \c d \ should be the data at the GUT scale. To 
evaluate the GUT scale data, RGE analysis of Yukawa couplings is necessary. This is the 
subject of the next section. 
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3 One-loop RGEs 



The data we need at the GUT scale are six quark masses, three angles and one CP-phase in 
the CKM matrix and three charged-lepton masses. In general, it is not a straightforward 
task to obtain the data at GUT scale by extrapolating the data at the electroweak scale. 
This is because RGEs of quark and charged-lepton Yukawa couplings include unknown 
neutrino Yukawa couplings. We can obtain a reliable answer only if neutrino Yukawa 
couplings are negligible. 

Fortunately, in models with the see-saw mechanism, the right-handed neutrinos are 
decoupled at the scale below the heavy Majorana masses, and the neutrino Yukawa cou- 
plings become irrelevant and disappear from RGEs of quark and charged-lepton Yukawa 
couplings. This fact makes our analysis straightforward, at least, below the heavy Majo- 
rana mass scale. In the following analysis, we use the RGEs independent of the neutrino 
Yukawa couplings at all scales between the electroweak and the GUT scales. Although 
this treatment fails at the scale above the right-handed Majorana neutrino masses, we 
expect that our results are reliable, since the heavy Majorana mass scale is close to the 
GUT scale. We will give a comment on this point in the last section. 

We analyze one-loop RGEs for gauge couplings and quark and charged-lepton Yukawa 
couplings in the MSSM. For simplicity, the universal soft SUSY breaking masses are 
assumed, msusY = Mz ~ 91GeV. In this case, the GUT scale is found to be Mq ~ 
2 x 10 16 GeV. The RGEs used in our analysis are as follows (HI- For gauge couplings, 



>h1, bi -J, (12) 



d(\lifl) 167T 

where hi = (3, —1, —33/5) for i — 3, 2, 1 corresponding to SU(3) C , SU F and U(l)y gauge 
couplings, respectively. We use the electroweak scale inputs such that (Mz) — sf/47T = 
0.118, a 2 (M z ) = # 2 /4tt = 0.0335 and a x {M z ) = g\j^ = 0.0168. On the other hand, 
RGEs for the Yukawa couplings of up- type quarks (Y u ), down-type quarks (Y d ), and 
charged-leptons (Y e ) are found to be 

( ' Sr 1 (P F S F + S F p F ), (13) 



<i(ln/i) IQtt 2 
where S F = Y F Y F with fermion index F = u,d, e, and 

A, = 35 u + ^+{tr[3SJ-(^ 3 2 + 3^ + ^ 1 2 )^ 1 



3x3 



3x3 



fa = 3S d + S u + {tr[3S d + S e )- (y^ 2 + 3^ + ^ 2 )|l 

Pe = ZS e +\tT[?>S d +S e } - (3^ + ^)}l 3 x3. (14) 
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To determine the input values of fermion Yukawa couplings at Mz, we have to fix tan/? 
by hand. For given tan/3, they are given by 



/ m 2 n , m 2 . m 2 \ / 9 „\ . . 

S u = diag^,^,^j(l+tan- 2 /3) , (15) 

S d = ^Mdiag^,^,^l4 M (l + tan 2 /3) , (16) 

S e = diag(^,^,^)(l + tan-^). (17) 

Here, we have worked out in the basis where S u and S e are diagonal. Given tan/3, it 
is easy to solve the above RGEs numerically. Although, as can be seen from the above 
RGEs, S u is no longer diagonal after running, we can extract the data we need at the 
GUT scale by the usual manner. For example, the CKM matrix at the GUT scale is given 
by Vkm = y^Ydi where V u and are the unitary matrices which diagonalize S u and Sd, 
respectively. 

Solving the GUT relation as discussed in the previous section by using the extrapolated 
data, we can obtain the light Majorana neutrino mass matrix, M u , as a function of a and 
cr through Eqs. ((SJ), (fTU|) and (fTTj) and the see-saw mechanism, M u = —M^M^Mc- Since 
this is the mass matrix at the GUT scale, it should be extrapolated to the electroweak 
scale according to the RGE for the effective dimension-five operator ^3] such that 



dM v 



l — {SJM„ + M u S e + (6tr [S u ] - 6g 2 2 - j^?) M„j . 



<i(ln/i) 16 

After the running, our resultant neutrino mass matrix are compared with the neutrino 
oscillation data. In the next section, the effect of the RGE running will be found to be 
small and almost negligible in our model. 



4 Numerical analysis and results 

Now we are ready to perform all the numerical analysis. Fist of all, note that the solution 
of the GUT relation exists, only if we take appropriate input parameters. This is because 
the number of our free parameters (fourteen) are almost the same as the number of inputs 
(thirteen). Therefore, in the following analysis, we vary two input parameters, m s and 
CP-phase 5 in the CKM matrix, within the experimental errors so as to find the solution. 
We take input fermion masses at Mz as follows (in GeV): 

m u = 0.00233 , m c = 0.677 , m t = 176 

m d = 0.00469 , m b = 3.00 

m e = 0.000487, m M = 0.103 , m T = 1.75 
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Here, the experimental values extrapolated from low energies to Mz were used |19j . For 
the CKM mixing angles in the "standard" parameterization, we input the center values 
measured by experiments as follows [20 : 



0.219 + 0.226 0.037 + 0.043 0.002 + 0.005 
s i2 — 5 ' 523 ~~ 2 ' Sl3 ~~ 2 ' 

Since it is very difficult to search all the possible parameter region systematically, we 
present only reasonable results we found. In the following, we show our analysis in detail 
in the case tan/5 = 45. Here we take m s = 0.0729(GeV) and 5 = 86.4°. In this case, the 
CKM matrix at Mz is given by 



V KM (M Z 



/ 0.975 0.222 0.000220 - 0.00349? 

-0.222 -0.000136i 0.974- 0.0000311? 0.0400 
V 0.00869 - 0.00340?' -0.0390 - 0.000777? 0.999 



After RGE evolutions to the GUT scale, fermion masses (up to sign) and the CKM matrix 
are found as follows (in GeV) : 

m u = 0.00103 , m c = 0.301 , m t = 134 
m d = 0.00171 , m s = 0.0265 , m b = 1.56 
m e = 0.000413 , m„ = 0.0872 , m T = 1.69 



and 



Vkm(M g ) 



( 0.975 0.222 0.000175 - 0.00279? 

-0.222 - 0.000121?' 0.974 + 0.000129? 0.0320 
V 0.00695 - 0.00272?' -0.0312 - 0.000626?' 0.999 



in the standard parameterization. 

These outputs at the GUT scale are used as input parameters in order to solve Eqs. (jHJ) 
and (jSJ). The contours of solutions of each equations are depicted in Fig. 1, where the 
signs of the input fermion masses have taken as — for m u , m c , md, and m s , and + for m t 
and rrib- We find a solution 

k = -0.0134 + 0.000786?' , 

\cd\ = 6.38 , (19) 

which leads to the reasonable results for the neutrino sector. 1 

Now we can describe all the fermion mass matrices as functions of a and cr by using 
the above solution through the mass matrix formulas of Eqs. (jSJ), (fTUj) and (fTTj) . Since 

1 Although there are other solutions, we do not address them. They lead to small sin 2 2#23 inconsistent 
with the atmospheric neutrino oscillation data, and are out of our interest. 
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the PMNS mixing angles are independent of cr, we can plot the angles as functions of a. 
The results are depicted in Fig. 2, and we can see that the resultant mixing angles are 
very sensitive to a. For a = 3.198[rad], we obtain sin 2 29u = 0.722, sin 2 2023 — 0.881 and 
sin 2 2^i3 = 0.164 at the GUT scale. After running this results back to the electroweak 
scale according to RGE of Eq. (fTK|) . we find 

sin 2 29 12 = 0.723, sin 2 26 23 = 0.895, sin 2 29 13 = 0.164 . (20) 

Note that RGE running effects are almost negligible. The ratio Am^/Am^ is also inde- 
pendent of cr, where Am 2 Q and Am^ are the oscillation parameters relevant for the solar 
and the atmospheric neutrino deficits, respectively. In Fig. 3, the ratio at the GUT scale 
is depicted as a function of a. After the RGE running, we find, at the electroweak scale, 

— -f = 0.188 . 21 
Am| 

RGE running effects are almost negligible also for the ratio. The neutrino mass matrix at 
the electroweak scale and the PMNS matrix which lead to the above results are as follows: 



M v = c R l 



I 14.7 + 48.0i -46.9 - 96.5i -446 - 370i 

-46.9 - 96.5i -347 + 90.3i 1030 + 702i ] (GeV) , 
\ -446 - 370i 1030 + 702i -497 - 914i 



and 



U 



PMNS 



( 0.168 + 0.838* -0.467 + 0.0940* -0.00508 + 0.207z 

0.0519 + 0.498i 0.651 - 0.0473i 0.0189 - 0.569i 
V 0.0745 + 0.116i 0.450 - 0.381i 0.431 + 0.669i 



They have been calculated in the basis where the charged-lepton mass matrix is diagonal 
and all the matrix elements are real and positive. The neutrino mass eigenvalues are 
determined by fixing cr appropriately. For example, taking cr = 3.01 x 10 13 leads to 
Amg = 3.76 x 10 _4 eV 2 and Am^ = 2.00 x 10~ 3 eV 2 . In this case, three mass eigenvalues of 
the right-handed neutrinos are found to be (in GeV) M Rl = 1.64 x 10 11 , M R2 = 2.50 x 10 12 
and Mr 3 = 8.22 x 10 12 , respectively. These results are in agreement with the results of 
the recent neutrino oscillation analysis |2I| with the large angle solution with Mikheev- 
Smirnov-Wolfenstein effect [22] for the solar neutrino deficit. 

Our results are very sensitive to the input values of m s and 5. As mentioned above, the 
solution of Eqs. (jHJ) and © exists only if the input values of m s and 5 are appropriately 
fixed. In Fig. 1, as m s is taken to be larger for fixed S, the contour of Eq.fjHJ) moves to the 
left, while the circle moves to the right. Thus, the solution disappears eventually. On the 
other hand, as 5 is taken to be larger for fixed m s , both of the contours move to the left, 
but the circle moves faster. Eventually, the solution disappears. This disappearance of 
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the solution occurs by several percent changes of the input parameters. Therefore, with 
fixed m s (5), 5 (m s ) is almost determined in order for the solution to exist. 

Next let us see a more important parameter dependence. When 5 and m s are fixed 
appropriately and the solutions are found, we obtain the similar plot as in Fig. 2 for the 
mixing angles in the PMNS matrix. However, the hight of the peaks depends on the input 
5 values. As in the past works [TT] [T3] , we mostly find that the resultant sin 2 26* 2 3 is too 
small to be consistent with the experimental data. Only the special values of the inputs 
and a can give the results being in agreement with all the neutrino oscillation data. The 
reasonable results we found are as follows: 



tan (3 


m s {M z ) 


5 


a 


sin 2 26 12 


sin 2 29 23 


sin 2 29 13 


Am^/Am^ 


40 


0.0718 


93.6° 


3.190 


0.738 


0.900 


0.163 


0.205 


45 


0.0729 


86.4° 


3.198 


0.723 


0.895 


0.164 


0.188 


50 


0.0747 


77.4° 


3.200 


0.683 


0.901 


0.164 


0.200 


55 


0.0800 


57.6° 


3.201 


0.638 


0.878 


0.152 


0.198 



The mass eigenvalues of the light Majorana neutrinos are determined by cr. In other 
words, cr is a function of the oscillation parameter, and appropriate value of cr is fixed 
so as to be suitable for the neutrino oscillation data. These relations are listed in the 
following. 



tan/3 


cr x m„ 3 (TeV) 


Cr X \ —fry 

K V io- 3 eV 


40 


2.18 


6.26 x 10 13 


45 


2.08 


6.02 x 10 13 


50 


1.81 


5.23 x 10 13 


55 


1.42 


4.04 x 10 13 



Here, m u:j is the mass eigenvalue of the heaviest light Majorana neutrino. As mentioned 
above, our resultant neutrino oscillation parameters are sensitive to all the input param- 
eters. In other words, if we use the neutrino oscillation data as the input parameters, 
the other input, for example, the CP-phase in the CKM matrix can be regarded as the 
prediction of our model. It is a very interesting observation that the CP-phases listed 
above are in the region consistent with experiments [2*U] . 

5 Other predictions 

Now all the mass matrices have been completely determined in our model as discussed 
in the previous section. There are other interesting physical observables which can be 
calculated by using the concrete fermion mass matrices in our model. 

The averaged neutrino mass relevant for the neutrino-less double beta decay [22] can 
be read off from (1,1) component of the light Majorana neutrino mass matrix M v . The 
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CP-violation in the lepton sector is characterized by the Jarlskog parameter [21] defined 

as 



Jcp = Im U e2 U; 2 U*JJ^ 



',22) 



where Ufi is the PMNS matrix element. 

It is well known that the SO(10) GUT model possesses a simple mechanism of baryo- 
genesis through the out-of-equilibrium decay of the right-handed neutrinos, namely, the 
leptogenesis [25J. The amount of the created baryon asymmetry is characterized by the 
CP-violating parameter e estimated as 



e ~ 



E 



Im 



(M D Ml 



2 

DJlj 



M Rl 



3=2,3 v 2 sin 2 



(3(M D Ml) n M Rj 



(23) 



Here, Md is the neutrino Dirac mass matrix in the basis where the right-handed Majorana 
mass matrix Mr is real and diagonal, and Mr. is the mass eigenvalue of the right-handed 
Majorana neutrino of the i-th generation. 

These quantities are evaluated by using the results presented in the previous section, 
and results are listed in the following. 



tan p 


(m u ) ee (eV) 


Jcp 


e 




40 


0.00122 


0.00110 


7.39 x 10" 


-5 


45 


0.00118 


-0.00429 


6.80 x 10" 


-5 


50 


0.00119 


-0.00631 


6.50 x lO- 


-5 


55 


0.00117 


-0.00612 


ll. 2 x 10- 


-5 



Here, cr was fixed so that Am^ = 2 x 10~ 3 eV 2 . Unfortunately, both of the averaged 
neutrino mass and the Jarlskog parameter may be too small to expect their evidences in 
future experiments. On the other hand, the CP-violating parameter is too large to be 
consistent with the observed baryon asymmetry. The problem of this too large baryon 
asymmetry can be easily avoided by assuming an inflationary universe whose reheating 
temperature is smaller than the right-handed neutrino masses. However, in this case, the 
leptogenesis scenario no longer works, and another scenario of baryogenesis such as the 
Affleck-Dine mechanism [2E] may be applicable. 

Sizable lepton-flavor violation (LFV) can be expected in SUSY models, and the LFV 
processes are one of the most important processes as the low-energy SUSY search |2Zj. In 
detailed analysis, concrete fermion mass matrices are necessary. It is worth investigating 
the LFV processes in our model. 



6 Summary and comments 

We have discussed the minimal SUSY SO(10) model, where only one 10 and one 126 
Higgs multiplets have Yukawa couplings with fermions. This model can determine all the 
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fermion mass matrices with only a few free parameters. It is known that, in the absence of 
CP violation, this model cannot incorporate the realistic neutrino mass matrix consistent 
with the neutrino oscillation data. We examined more general case in which CP-violating 
phases in the fermion mass matrices were introduced. In this case, using the experimental 
data of six quark masses, three angles and one CP-phase in the CKM matrix and three 
charged-lepton masses, the light Majorana neutrino mass matrix can be determined as 
the function of only one free parameter a through the GUT relation and the see-saw 
mechanism. Here we do not count cr, which plays only the role to fix the overall scale of 
the neutrino masses. In order to connect the mass matrix data between the electroweak 
and the GUT scales, we have analyzed the one-loop RGEs of the charged fermion Yukawa 
couplings and the effective dimension-five operator of the light Majorana neutrinos. 

We found that there was the parameter region in which the predicted neutrino mass 
matrix can be consistent with the current neutrino oscillation data. This parameter region 
is severely constrained, and all the parameters, even the CP-phase in the CKM matrix, are 
almost fixed for given tan/3. In other words, we can regard the neutrino oscillation data 
as the inputs while the CP-phase as the output. Interestingly, our results consistent with 
the neutrino oscillation data are also consistent with experimental results of CP-violation 
in the quark sector. 

In our RGE analysis, we have used RGEs introduced in Sec. 3 at all the energy scales 
between the electroweak scale and the GUT scale. This treatment fails at the energy larger 
than the right-handed Majorana neutrino mass scale, where the right-handed neutrinos 
are no longer decoupled. Above that scale, we have to take RGE of neutrino Yukawa 
coupling into account. However, as found in Sec. 4, the Majorana mass scale is about 
10 13 GeV, and close to the GUT scale in the sense of the RGE (logarithmic) running. 
Thus, we expect that our results can be still reliable. In fact, we can find that changes 
of our input values at the GUT scale remain within several percents even if the neutrino 
Yukawa coupling RGE is taken into account, and that the results presented in Sec. 4 are 
almost unchanged. 

In addition, we have neglected the SUSY threshold corrections in our RGE analysis. 
It is known that with large tan(3 the down-type quark mass matrix is potentially affected 
by large SUSY threshold corrections and we can neglect them with a limited region of 
the soft SUSY breaking parameters [2H]- If the SUSY threshold corrections are large, 
it should be taken into account in our analysis, and will change the input values at the 
GUT scale and eventually the final result. It is an interesting question whether the result 
gives a better or worse fit, and it is worth analyzing the RGEs taking SUSY threshold 
corrections into account. 

When the CHOOZ data [2H| is concerned, our results seems to be in sever situation 
[3*U] . However, more comprehensive analysis [21] gives rather loose constraint. In any case, 
our results lie in the allowed region at 99% C.L. in the global analysis [2*T] . and we have 
concluded that they are in agreement with the current neutrino oscillation data. 2 This 

2 We thank Osamu Yasuda for his analysis concluding that our results is 2.6 a away from the current 
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situation seems to be inevitable in our model, since two angles, sin 2 2# 2 3 and sin 2 2# 13 , 
are strongly correlated and have the peaks at almost the same a value as can be seen in 
Fig. 2. This can be easily avoided, when we extend the minimal model and introduce new 
Higgs multiplets. An extended model has new free parameters in fermion mass matrices, 
and can be in much better agreement with the neutrino oscillation data, larger sin 2 2#23 
but smaller sin 2 2# 13 . However, what we have found in this paper is that the minimal 
SO (10) model, the model with only minimal set of Higgs multiplets, is still viable without 
any extention. 

We have concentrated our discussion only on the Yukawa sector with one 10 and 
one 126 Higgs multiplets. In order to complete our model, it is necessary to construct 
a concrete Higgs sector which can realize all the assumptions in our model, the correct 
GUT symmetry breaking, doublet-doublet (triplet) Higgs mass splitting etc. Such a Higgs 
sector may include new Higgs multiplets which affect the fermion mass matrices. If it is the 
case, our model becomes less predictive, although the model can be in better agreement 
with the neutrino oscillation data. 3 This issue is highly non-trivial, and a further work is 
needed. 

Although, in this paper, the SUSY SO(10) model has been considered, non-SUSY 
SO(10) model with intermediate symmetry breaking scales is also worth investigating. In 
this case, the strategy is completely analogous to that we discussed in this paper, except 
that RGEs are replaced by the ones of non-SUSY case. 
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Figure 1: Contour plot on complex K-plane. The vertical line and the circle correspond 
to the solutions of Eqs. (jHJ) and (jHJ), respectively. 
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Figure 2: Three mixing angles in the PMNS matrix as functions of er[rad]. The graphs 
with the highest, middle and lowest peaks correspond to sin 2 26* 2 3, sin 2 20 12 and sin 2 2# 13 , 
respectively. The plots of sin 2 2^23 arid sin 2 26*i3 have the sharp peaks at o ~ 3.2 [rad], 
while sin 2 2#i2 has the sharp peak at a ~ 3.3 [rad]. 
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Figure 3: The ratio (Log 10 [AmQ/ Am^]) as the function of cx[rad]. 
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